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^ ' Abstract 



Levy copulas are the most general concept to capture jump dependence in 
multivariate Levy processes. They translate the intuition and many features 
of the copula concept into a time series setting. A challenge faced by both, 
distributional and Levy copulas, is to find flexible but still applicable mod- 
els for higher dimensions. To overcome this problem, the concept of pair 
copula constructions has been successfully applied to distributional copulas. 
In this paper, we develop the pair construction for Levy copulas (PLCC). 
Similar to pair constructions of distributional copulas, the pair construction 
of a (i-dimensional Levy copula consists of d(d — l)/2 bivariate dependence 
functions. We show that only d — 1 of these bivariate functions are Levy 
copulas, whereas the remaining functions are distributional copulas. Since 
there are no restrictions concerning the choice of the copulas, the proposed 
pair construction adds the desired flexibility to Levy copula models. We dis- 
cuss estimation and simulation in detail and apply the pair construction in 
a simulation study. 
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1. Introduction 

Many financial and nonfinancial applications need multivariate models 
with jumps where the dependence of the jumps is captured adequately. To 
this end, Levy processes have been applied in the literature. However, al- 
though the recently introduced concept of Levy copulas enables modeling 
the dependence in Levy processes in a multivariate setup, known parametric 
Levy copulas are very inflexible in higher dimensions, i.e., they consist of 
v ery few pa—. In this paper, we show that, shnUar to the pah copuia 
construction of distributional copulas going back to Joe |13| . Levy copulas 
may be constructed from a constellation of parametric bivariate dependence 
functions. Because these dependence functions may be chosen arbitrarily, 
the resulting Levy copulas flexibly capture various dependence structures. 

Levy processes are stochastic processes with independent increments. 
They consist of a Brownian motion part and jumps. Due to the jumps, Levy 
processes capture stylized facts observed in financial data as non-normality, 



14]). At the same time, 



excessive skewness, and kurtosis (see, e.g., Johannes 
they stay mathematically tractable and allow for derivative pricing by change 
of measure theory. For these reasons, intensive research is conducted on the 
statistical inference of Levy processes (see, e.g., Lee and Hannig [18] and the 
references therein). 

The fundamental work for multivariate applications of Levy processes is 



the seminal paper of Kallsen and Tankov 



16] , where the concept of Levy 



copulas is introduced. This concept transfers the idea of distributional cop- 



ulas to the context of Levy processes. Distributional copulas (normally just 
referred to as copulas) are functions which connect the marginal distribu- 
tion functions of random variables to their joint distribution function. They 
contain the entire dependence information of the random variables (see, e.g., 



Nelsen 19| for an introduction to copulas). In the same sense, the theory of 
Levy copulas enables to model multivariate Levy processes by their marginal 
Levy processes and to choose a suitable Levy copula for the dependence 
structure separately. For papers regarding the estimation of Levy copulas in 
multivariate Levy processes and app lications see, e.g., the recent papers of 
Esmaeili and Kliippelberg Jj], [nj, [llj and references therein. 

All papers involving Levy copulas focus on rather small dimensions since 
higher-dimensional flexible Levy copulas are difficult to construct. A similar 
effect has been observed during the first years of literature on distributional 
copulas, where mainly 2-dimensional distributional copulas have been ana- 
lyzed. One solution regarding distributional copulas has been the develop- 



ment of very flexible pair constructions of copulas going back to Joe 13 1 
and further developed in a series of papers (see, e.g., Bedford and Cooke |5| 
or Aas et al. jlj|). In pair copula constructions, a <i-dimensional copula is 
constructed from d(d — l)/2 bivariate copulas. Here, d — 1 of the bivariate 
copulas model the dependence of bivariate margins, whereas the remaining 
bivariate copulas model certain conditional distributions, such that the entire 
(i-dimensional dependence structure is specified. 

Levy copulas are conceptually different from distributional copulas. While 
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^-dimensional distributional copulas are distribution functions on a [0, l] d hy- 
percube, d- dimensional Levy copulas are defined on IR d and relate to Radon 
measures. Therefore, the idea of pair constructions for copulas is not directly 
transferable to Levy copulas and up to now it has not been clear whether it 
is possible at all. In this paper, we show that a pair copula construction of 
Levy copulas (PLCC) is indeed possible. It also consists of d(d — l)/2 bivari- 
ate dependence functions but only d — 1 of them are Levy copulas, while the 
remaining ones are distributional copulas. For statistical inference, we derive 
sequential maximum likelihood estimators for an arbitrary pair construction 
of Levy copulas as well as a simulation algorithm. We analyze the applica- 
bility of the concept in a simulation study. The estimation and simulation 
algorithms show encouraging results in a finite sample setting. 

The remainder of the paper is structured as follows. In Section |2~T| we re- 
view the theory of copulas for random variables and pair copula constructions 
of such copulas. In Section 12. 2[ we address the theory of Levy processes and 
Levy copulas. Our pair construction of Levy copulas is derived in Section [3j 
In Section HI we provide simulation as well as maximum likelihood estimation 
methods. Section [5] contains simulation studies probing the simulation and 
estimation algorithms in finite samples and Section [6] concludes. 

2. Preliminaries 

In this section, we briefly recall necessary theory on copulas, pair copulas, 
Levy processes, and the Levy copula concept. 
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2.1. Copulas and Pair Copula Construction 

Let X = (X\, . . . , X d ) be a random vector with joint distribution function 
F and continuous marginal distribution functions Fi, i = 1, . . . , d. The copula 
C of X is the uniquely defined distribution function with domain [0, l] d and 
uniformly distributed margins satisfying 

F(xt, ...,x d ) = C[Fi(xi), . . .,F d (x d )]. 

By coupling the marginal distribution functions to the joint one, the copula 
C entirely determines the dependence of the random variables X\, . . . ,X d . 
While many 2-dimensional parametric families of copulas exist, see, e.g., 
Nelsen 19|, the families for the d- dimensional case suffer from lack of flex- 



ibility. To overcome this problem, the concept of pair copula construction 

has been developed (see, e.g., Joe jl3| for the seminal work or the detailed 

I 

introductions in Aas et al. Ilk Bedford and Cooke 51, and Berg and Aas 

n 

|7|). In a pair copula construction, a <i-dimensional copula C(u\, . . . ,u d ) is 
constructed of d(d — l)/2 bivariate copulas. Of these bivariate copulas, d— 1 
bivariate copulas directly model d — 1 2-dimensional margins of the copula 
C, whereas the other bivariate copulas indirectly specify the remaining parts 
in terms of conditional distributions. Since the number of possible combina- 
tions grows rapidly with the dimension, Bedford and Cooke ||], 0| introduced 
a graphical model, called regular vines (R- vines), to describe the structures 
of pair copula constructions. 
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T2 (12) ( 23 

Cl,3|2 

Figure 1: Example of a pair copula construction of a 3-dimensional copula. It consists of 
2 trees and 3(3-l)/2=3 bivariate copulas Cx t 2, C23 and C\^\2- 



An example of a regular vine for the 3-dimensional case is given in Figure 
[TJ It shows three dimensions (labeled 1,2 and 3) and two trees (labeled Tl and 
T2) of dependence functions. The first tree (Tl) contains the two bivariate 
copulas, C%,2 and 6*2,3, modeling the dependence between dimensions 1 and 
2 and dimensions 2 and 3, respectively. Thus, tree Tl completely determines 
these two bivariate dependence structures. It also indirectly determines parts 
of the dependence between dimensions 1 and 3, but not necessarily the en- 
tire dependence. For instance, if the pairs 1,2 and 2,3 are each correlated 
with 0.9, then 1 and 3 cannot be independent but their exact dependence 
is not specified. In particular, the conditional dependence of 1 and 3 given 
2 is not specified. Therefore, in the second tree (T2), this bivariate con- 
ditional dependence is modeled with another copula, Cx^\2- Together, the 
three bivariate copulas fully specify the dependence of the three dimensions. 
Since the choice of all three bivariate copulas is arbitrary, the vine structure 
provides a very flexible way to construct multidimensional copulas. In the 
rf-dimensional case, d(d — l)/2 bivariate copulas are needed and arranged in 



d — 1 trees (see, e.g., Joe 



13]). There are special cases of regular vines, e.g., 
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C- vines or D- vines (see, e.g., Aas et al. [l| for a more detailed introduction). 

2.2. Levy Processes and Levy Copulas 

Detailed information about Levy processes may be found in Rosinski 



20], Kallenberg [15] or Sato 



21] . Introductions to Levy copulas are given in 
Kallsen and Tankov [16| or Cont and Tankov [8]. Here, we give a very short 
overview of both. 

Let (O, P) be a probability space. A Levy process (L t ) tg R + is a 
stochastic process with stationary, independent increments starting at zero. 
Levy processes can be decomposed into a deterministic drift function, a Brow- 
nian motion part and a pure jump process with a possibly infinite number 



of small jumps, see, e.g., Kallenberg 



15], Theorem 15.4 (Levy Ito decompo- 



sition). In this paper, we focus on spectrally positive Levy processes, which 
are Levy processes with positive jumps only. This facilitates the notation 
considerably and in many relevant cases it is sufficient to consider positive 
jumps only. However, all results of the paper may be extended to the general 
case. The characteristic function of the distribution of such an Revalued 
spectrally positive Levy process L+, at time t, is given by the Levy-Khinchin 
representation (see Kallenberg [15] ) 

r l, { ~ ) ^{i ( /(-,.:> -~(z,Zz)+ I (e i{z ' x) - l)"(dx) I }. (1) 



expjt ^i(j,z) - ~(z,Ez) + j 



Here, 7 G R corresponds to the drift part of the process and E is the 
covariance matrix of the Brownian motion part at time t — 1. The Levy 



measure v is a measure on M. d which is concentrated on the positive domain 
\ {0} with j Rd xu(dx) < oo. The Levy measure completely characterizes 
the jump parts of the Levy process, where v(A) for A G S(Rl) is the expected 
number of jumps per unit of time with jump sizes in A. A spectrally positive 
Levy process with positive entrees of 7 and £ = is called subordinator. It 
has no negative increments. 

An interesting example for a one-dimensional subordinator is the stable 
subordinator. It is heavy tailed and therefore suggested as a loss process for 



operational risk models. In Basawa and Brockwell 
a stable subordinator on IR + is defined by 



3|, the Levy measure of 



u(B) 



l^U) — -az, 



where a G (0, 1) and > 0. 

Related to the Levy measure, its tail integral is defined by (see, e.g., 
Definition 3.1 in Esmaeili and Kliippelberg [9]) 



U(X!, ...,x d )= < 



f([xi, 00) X 




00 



x [x d , 00)) if (xi, . . . , x d ) G [0, oo) d \{0}, 
if Xi = 00 for at least one i, 
if (xi, ...,x d )=0. 



The tail integral U of a spectrally positive Levy process uniquely deter- 
mines its Levy measure v. We define the marginal tail integrals Uk for any 
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dimension k = l,...,d of the multivariate Levy process in a similar way. 
For one- dimensional spectrally positive Levy measures u, the tail integral is 
U(x) = z/([x, oo)), i.e., the expected number of jumps per unit of time with 
jump sizes larger or equal to x. For the one-dimensional stable subordinator, 
the tail integral can be explicitly calculated and inverted for x, u > 0, 



The inverse of the tail integral is needed for the simulation of the process. 

Dependence of jumps of a multivariate Levy process can be described 
by a Levy copula which couples the marginal tail integrals to the joint one. 
A o?-dimensional Levy copula is a measure defining function €(ui, . . . ,Ud) : 
[0, oo] d —> [0, oo] with margins <£k(uk) '■= £(oo, . . . , oo, Uk, oo, . . . , oo) = Uk 
for all Uk € [0, oo] and k — 1, . . . , d. In particular, let U denote the tail integral 
of a spectrally positive <i-dimensional Levy process whose components have 
the tail integrals U\, . . . , Ud- Then, there exists a Levy copula € such that 
for all (xi, . . . , Xd) € M^j_ 



Conversely, if £ is a Levy copula and U\, . . . ,Ud are marginal tail integrals 
of spectrally positive Levy processes, Equation (j2J) defines the tail integral 
of a (i-dimensional spectrally positive Levy process and Ui, . . . , Ud are the 
tail integrals of its components. Both statements are often called the Sklar's 




U(xi,...,x d ) = ^{U^xx), . . . ,U d {x d )). 



(2) 
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theorem for Levy copulas and are proved, e.g, in Cont and Tankov 8|. 

In this paper, we focus on Levy copulas for which the following assump- 
tion holds. 

Assumption 1. 

Let Cl,...^ be a Levy copula such that for every I C {1, . . . , d} nonempty, 



lim <£i,...,d(ui, ... ,u d ) = £i,...,d(wi, ••• ,u d ) \( Ui ) ieI =oo- (3) 

[Ui)i &I -*(X> 



This is a rather weak assumption on the Levy copula and is assumed in 



many papers, e.g., in Tankov 



22fl . It means that the Levy copula has no new 



information at the points Ui = oo which is not already contained in the limit 
for Ui — > oo. We need it since it ensures a bijection between a Levy copula on 
ffifj. and a positive measure /xi,...^ on with one-dimensional Lebesgue 

margins. This measure is given by 

A*w((M]) = Ve a d (MD, (4) 

where a, b e with a < b, component-wise, and d refers to the £i v .. i( r 
volume of the d-box [a, b] which is defined as 

Vc 1 ,..J[a,b]) = J2^(c)^i,...Ac)- 
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The sum is taken over all vertices c of [a, b] and 



sgn(c) 



1 if Ck = afc for an even number of k, 
1 if Ck = ak for an odd number of k. 



Furthermore, any positive measure //i,...^ on M~L with Lebesgue margins 
uniquely defines a Levy copula on that satisfies Assumption [1] by 

• • • , u d ) := «l,...,d([0, Ml] x ... x [0, u d }) 

and by setting 

£i,...,d(ui, u d ) \{ Ui ) ieI=00 := lim /ii,...,d([0, Mi] x ... x [0, u d \). 



131- 



These results are proved, e.g., in Section 4.5 in Kingman and Taylor 

An example for a Levy copula which is used later in the paper is the 
Clayton Levy copula. For spectrally positive, 2-dimensional Levy processes 
it is given on M + by 

£(u,v)= ( u - s + v- y 1/9 . (5) 

Here, 6 > determines the dependence of the jump sizes, where larger values 
of 9 indicate a stronger dependence. 
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3. Pair Levy Copulas 

In this section, we present the pair construction of d- dimensional Levy 
copulas. In particular, we show that analogously to the pair construction of 
distributional copulas, d(d — l)/2 functions of bivariate dependence may be 
arranged such that they define a d- dimensional Levy copula. In Sections 13.21 
and 13. 'S\ we provide illustrating examples how to construct multivariate pair 
Levy copula constructions. Readers not interested in the technical parts may 
read these examples first. 

3. 1 . Technical Part 

The central theorem for the construction is Theorem [2j It states that two 
(d — 1) -dimensional Levy copulas with overlapping (d — 2)-dimensional mar- 
gins may be coupled to an <i-dimensional Levy copula by a new, 2-dimensional 
distributional copula. Ensured by vine constructions (see Bedford and Cooke 
6|) and starting at (d— 1) = 2, Theorem [2] therefore enables to sequentially 
construct Levy copulas out of 2-dimensional dependence functions, i.e., 2- 
dimensional distributional copulas and Levy copulas. Before we state the 
theorem, for convenience, we recall some definitions which can be found, 
e.g., in Ambrosio et al. |2|. 

Definition 1. 

A positive measure on (IR+, i3(R+)) that is finite on compact sets is called a 
positive Radon measure. 

Let (X, S) and (Y, J 7 ) be measure spaces and let f : X — > Y be a measurable 
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function. For any measure /i on (X, £), we define the Push Forward Measure 
f#fi in (Y, J 7 ) by 

:= /x (r\K)) \/K G J 7 . 

Let fx be a positive Radon measure on and x i— > £ x a function which 
assigns a finite Radon measure £ x on to each x G Mi . We say this map 
is yU-measurable if x H- £ X (B) is fi-measurable for any B G £>(M™). 

Definition 2 (Generalized Product). 

Let n be a positive Radon measure onM.+ andx i-> £ x a ^-measurable function 
which assigns a probability measure £ x on to each x G Mi.. We denote by 
[L (g) £ x the Radon measure on M+ +m defined by 

H®U B )-= [ if tB(x,y)dUy) ) dfi(x) WBeB(KxR^), 

where K C is any compact set. 

We also need a theorem which states that a Radon measure may be 
decomposed into a a projection onto some of its dimensions and a probability 
measure. For a proof see Theorem 2.28 in Ambrosio et al. 
sentence after Corollary 2.29 there. 



2| and also the 



Theorem 1 (Disintegration). 

Let /ii^.^+m be a Radon measure on IR+ +m , 7r : M+ +m i— y the projection 
on the first d variables and //i,... )t j = n^fii^.^d+m- Let us assume that /ii,...^ is 
a positive Radon measure, i.e., that /ii i ... jC ; +m (i^ x M.™) < oo for any compact 
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set K C R+. Then, there exists a finite measure £ x in R™ such that a;4^ 
is /ii v „ ^-'measurable, £ x is a probability measure almost everywhere in 
and 

/ l B (x,y)dn h ... !d+m (x,y) = / / l B (x,y)£ x (y) \ d/i li ... 4 (x), 

this is pL 1: _ :d+m (B) = Hi,...,d®£x{B) for any B E B(K x R™), w/iere i^Clf 
is any compact set. 

We are now able to state the main theorem. 
Theorem 2 (Pair Levy Copula Composition). 

Let and £2,...^ &e two Levy copulas on R+ 1 where €i r ,, t d-i is a Levy 

copula on the variables Ui, . . . ,u d -i and (£2,...^ ^ s a ^ewy copula on the vari- 
ables «2, • • • , Md. Denote the corresponding measures on R+ _1 6y yLii,...,d— 1 
H2,...,d, respectively. Suppose that the two measures have an identical (d — 2)- 
dimensional margin /i2,...,d-i on the variables u 2 , . . . , Ud-i- Then, we can 
define a Levy copula on R^ by 

ei,..,d(ui,...,Ud) := J ^(Fxi^^.^^Mi),^!^,...^^^^)^,...^-!^,.-.,^-!), 

[0,u 2 ]x...x[0,u d _i] 

where Fi^...^^ is the one- dimensional distribution function corresponding 
to the probability measure £i| U2) ..., Ud _i /ram £/ie decomposition of Hi t ...,d-i into 

fJ>l,...,d-l = fJ>2,...,d-l <8> ^l|u 2 ,...,W d _i) 
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Fd\u 2 ,....u d _ 1 is the one- dimensional distribution function corresponding to the 
probability measure £d\u 2 ,-,u d _i from the decomposition of ^2,...,d into 

^2,...,d = ^2,...,d-l <8> £d|ua,— >«d-l' 

and C is a distributional copula. Since Levy copulas are functions on we 
set for every I G {1, . . . , d} nonempty, 

■ ■ ■ , «d)lK) ieI =oo := , lim ■ ■ ■ , u d ). (6) 

The theorem, which is proved in the appendix, illustrates how to construct 
a ci-dimensional Levy copula from two (d— l)-dimensional Levy copulas with a 
common margin. Applying the theorem recursively, these (d— l)-dimensional 
Levy copulas can be constructed from (d — 2)-dimensional ones. This can be 
repeated down to construct 3-dimensional Levy copulas from bivariate ones. 
In higher dimensions, there are many ways for this procedure due to possible 
permutations of the dimensions and numerous possible pairwise combina- 
tions within the trees. The graphical visualization of the different resulting 
structures of pair construction is possibly by the concept of regular vines as 
developed in Bedford and Cooke |6|. Regular vines also help to construct 
pair Levy copulas top-down. This means to start with d — 1 bivariate Levy 
copulas and to combine them successively to 3,4,5,..., rf-dimensional Levy 
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Tl 





® 



T2 




® 



Figure 2: Pair construction of a 3-dimensional Levy copula out of 3(3 — l)/2 = 3 bivariate 
dependence functions. The functions £^2 and £2,3 in the first tree are Levy copulas, while 
Ci,3|2 i n the second tree is a distributional copula. 

copulas. The regular vine approach ensures that at each step the involved 
Levy copulas have sufficiently overlapping margins, and that therefore the 
theorem can be applied. To illustrate this procedure, we give two detailed 
examples. The first example refers to the most simple 3-dimensional 
Levy copula. The second, 4-dimensional example then illustrates how to 
sequentially add dimensions to the pair copula construction. 

3.2. Example: 3-dimensional Pair Levy Copula Construction 

A 3-dimensional example can be constructed applying Theorem [2] to com- 
bine two 2-dimensional Levy copulas by a distributional copula. As in the 
usual pair copula construction for distributional copulas, in Figure [2] we use 
the vine concept to visualize the resulting dependence structure. The bi- 
variate dependence structures in the first tree are Levy copulas, whereas the 
copula in the second tree is a distributional copula. From Theorem [2] follows 
that 





[0,u 2 ] 
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is a Levy copula, where F 1 | U2 (m 1 ) is the one- dimensional distribution function 
corresponding to the probability measure £i| U2 from the decomposition of //i i2 
into 

^1,2 = f*2 ® 6|« 2 ( 7 ) 

and -F31U2 is the one-dimensional distribution function corresponding to the 
probability measure £ 3 | U2 from the decomposition of yU 2 ,3 into 

A*2,3 = A < 2®6|u a - (8) 

Remember that yU 12 is the Radon measure corresponding to £12. With 
Theorem [1] and the considerations after Assumption [1] we see that /i 2 in 
Equation (J7]) is the Lebesgue measure. Analogously, /X2,3 is the Radon mea- 
sure corresponding to (£ 2 ,3 and therefore /x 2 in Equation (jSJ) is the Lebesgue 
measure as well. 
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To check whether £1,2,3(^1,^2,^3) has the desired margins, we calculate 



£1,2,3(^1,^2, 00) 



Ci,3\ 2 (F 1 \ Z2 (u l ),F 3lz2 (oo))dz 2 
[0,112] 

= J C lt3 \ 2 (F 1 \ Z2 (u 1 ),l)dz 2 

[0,u 2 ] 

= J F llz2 (u 1 )dz 2 



[0,112] 



( 



[0,u 2 ] \[0,«i 



dz 2 



^1,2(21,2:2) 



[0,ui] x [0,u 2 ] 
£1,2(^1,^2). 



A similar procedure shows that 



£1,2,3(00,^2,^3) = £2,3(^2, u 3 ). 



As expected, we do not get such a direct representation of the third bivariate 
margin 

£1,2,3(^1,00,^3)= J Ci^ 2 (Fi\ Z2 {u 1 ),F 3 \ Z2 {u z ))dz 2 
[0,00) 

because this margin is not only influenced by the distributional copula C\^\ 2 
but also by £ 12 and £2,3. However, we can adjust the bivariate margin of the 
first and third dimension by changing C\ 312 without affecting the other two 
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T1 © z © z © 

^2,3 ^3,4 



T2 (23) (34 

C*2,4|3 

Figure 3: Pair construction of the second three dimensions of a 4-dimensional Levy copula 
out of 3(3 — l)/2 = 3 bivariate dependence functions. The functions £2.3 and £3.4 in 
the first tree are Levy copulas, while C 2j 4|3 in the second tree is a distributional copula. 
The Levy copula £2,3 is the same Levy copula as in Figure [5] which refers to a the pair 
construction of the first three dimensions. 

T1 © ; © 1 © ; 

^1,2 ^2,3 ^3,4 



T2 (12) (23J (34 

C*l,3|2 C2,4|3 

T3 Cl§D (243) 

Cl,4|2,3 

Figure 4: Combination of the first three dimensions and the second three dimensions to a 
pair construction of a 4-dimensional Levy copula. It consists of 4(4 — l)/2 = 6 bivariate 
dependence functions. Only the functions in the first tree are Levy copulas, while the 
functions in the second and third tree are distributional copulas. 



bivariate margins. 

3.3. Example: 4-dimensional Pair Levy Copula Construction 

Considering 4 dimensions, we need two 3-dimensional Levy copulas with 
an identical 2-dimensional margin. Here, we reuse the Levy copula from 
Example 13.21 for the first three dimensions. The second 3-dimensional Levy 
copula is constructed in the same way and has the vine representation shown 
in Figure [3j 

Notice that the Levy copula £2,3 is used in both 3-dimensional pair Levy 
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copulas. Therefore, the marginal Levy copulas 

£1,2,3(00,^2, = £2,3(^2,^3) = £2,3,4(^2,^3, 00) 

are the same and we can apply Theorem [2] to construct a 4-dimensional Levy 
copula with the vine representation shown in Figure H] and 

£1,2,3,4(^1, u 2 ,u 3 ,U4,) = J Ci A \2, 3 (Fi\z 2 ,z a (ui), F 4 \ Z2:Z3 (u 4 ))dfi 2 ,3(z2, Z3) 

[0,«a]x[0,u 3 ] 

where Fx\ U2jU3 is the one-dimensional distribution function corresponding to 
the probability measure £i| U2 ,u 3 from the decomposition of 7x1,2,3 from the first 
pair Levy copula £1,2,3 into 

7^1,2,3 — 7^2,3 <S> £l|U2,«3' 

The one- dimensional distribution function i*4i U2 , U3 corresponds to the proba- 
bility measure ^4| U2 , U3 from the decomposition of 7x2,3,4 from the second pair 
Levy copula £2,3,4 into 

7*2,3,4 = 7^2,3 ® ^4|u 2 ,M3- 

4. Simulation and Estimation 

In this section we discuss the simulation of multivariate Levy processes 
as well as the maximum likelihood estimation of the pair Levy copula. We 
need the following assumption which is fulfilled by the common parametric 
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families of the bivariate (Levy) copulas. 
Assumption 2. 

In the following, we assume that all bivariate distributional and Levy copulas 
are continuously differentiable. 

4-1. Simulation 

The simulation of multivariate Levy processes built upon Levy copulas 
bases on a series representation for Levy processes and the following theorem. 

Theorem 3. 

Let v be a Levy measure on satisfying f R d(\\x\\ A l)dv{x) < 00, with 
marginal tail integrals Ui, i — 1, . . . , d, Levy copula (ti,...^ with corresponding 
measure /ii,...,d. Let (Vi)^ be a sequence of independent and uniformly [0, 1] 
distributed random variables and (Tj, . . . , r^ -1 ) ieN be a Poisson point process 
on with intensity measure >Lti,...,d_i from the decomposition of 

^i,...,d = fj>i,...,d-i <s> £d\ ui,..., u d _i, 

with ^d| tt i,...,tt d _i being a probability measure. For any value of Fj, . . . , F^ 1 , 
we suppose that Ff is a random variable with probability measure £ d | r i r d-i . 
Then, the process (L], . . . , L d ) te [ 0;1 ] defined by 

00 

^E^H^mM, J = l,...,d 

i=l 
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is a d- dimensional Levy process (L t ) f6 r 0> i] without a Brownian component and 
drift. The Levy measure of L t is v. 



Proof: The proof is similar to the proof of Tankov 



22j, Theorem 4.3. 



In practical simulations, the sum cannot be evaluated up to infinity and 
one omits very small jumps. The sequence (T^)^ is therefore only simulated 
up to a sufficiently large N, resulting in a large value of (see Rosinski 2oj] 
for this approximation). Note that large values of correspond to small 
values of the jumps (Tj), since the tail integral is decreasing. 

Based on the pair copula construction of the Levy copula, Tf, . . . ,Tf can 
be drawn conditionally on Tj in a sequential way. For convenience, assume 
that the pair Levy copula has a D-vine structure and that the dimensions 
are ordered from left to right. The dependence between Fj and T*j is then 
determined in the first tree of the pair construction by the bivariate Levy 
copula £1,2, and the distribution function F 2 | r i of Tj given Tj is derived in 
the following Proposition. 

Proposition 1. 

Let £ 12 be a 2-dimensional Levy copula with corresponding measure fi 12 - 
Then, we can decompose 

Ati,2 = A*i <8> 6|«D 

where ^2\m is a probability measure and the distribution function for almost 
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all U\ G [0, oo) is given by 

F 2 \ U1 {U2) = • 



Proof: This is a special case of Tankov 22], Lemma 4.2. 

Inverting this distribution function allows the simulation of T^. Now sup- 
pose that we have already simulated the variables T 1 , . . . , d > 3 and 
we want to simulate the last variable T d . We already know from Theorem [1] 
that the distribution of the last variable, given the first d — 1, is a specific 
probability distribution and therefore we are interested in the correspond- 
ing distribution function F ( i| Ul) _ j „ d _ 1 . Having found Fd\ Ult ,„ !Ud -i, we can again 
invert it and easily simulate a realization of a random variable with this dis- 
tribution function. The next proposition provides i*di Ul within the pair 
construction of the Levy copula. 

Proposition 2. 

Let d > 3 and <£i d be a pair Levy copula, fii d the corresponding measure, 
7T the projection on the first d — 1 variables, and fix d-i — ^^/^l,...^ the push 
forward measure. Then, we can decompose 

V>X,...,d = A*l,...,rf-1 <E> £d\ui,...,u d -i, 
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where £d\ui,...,ud-i ^ s a probability measure on M + with distribution function 

9Cl,d|2,...,d-l(-Pl|ii 2 ,-,'Ud-i( U l)> Fd\u 2 ,...,u d _ 1 {Ud)) 



Fd\uu-,u d -i\ U d) 



dF 1 \ U2) _ Ud _ 1 {u l ) 



it 



l^i,...,d-i- almost everywhere. Moreover, i^i Wli .„ )U ._ 1 is continuously differen 
tiable. 

The proposition is proved in the appendix. Similar to Aas et al. l| 
shows how we can iteratively evaluate and invert the distribution function 

Fd\u 1 ,...,u d _ 1 - 

4-2. Maximum Likelihood Estimation 

It is usually not possible to track Levy processes in continuous time. 
Therefore, we have to choose a more realistic observation scheme. In the 
context of inference for pure jump Levy processes, a common assumption 
is that it is possible to observe all jumps of the processes larger than a 
given e (see, e.g., Basawa and Brockwell (1978, 1980) 3|, |4| or Esmaeili and 
Kliippelberg [9|]). 



ll| . we estimate the marginal Levy 



Following Esmaeili and Kliippelberg 
processes separately from the dependence structure. That is, we use all 
observations with jumps larger than e in a certain dimension and estimate 
the parameters of the one-dimensional Levy process. 

For the estimation of the dependence structure, i.e., the Levy copula, we 
can use the fact that the process consisting of all jumps larger than e in all 
dimensions is a compound Poisson process. We suppose that all densities 
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fx , ■ ■ ■ , fd of the marginal Levy measures exist and we denote the parameter 
vectors of the Levy copula and the marginal Levy measures by 5, 71, ... , 74, 
respectively. The likelihood function is given by 



L £ {S,jx,...,j d ) = 

Ml, 

e " Al '"" d * I f LM^i' 7i) • • • • • fd{xid, ld)Cx,...,d(Ux(xix, 71), . . . , Z7d(^id, 7d)> 



8=1 



where d = £x,...,d(Ux(e, 71), . . . , Ud(s, 7^), S) , iV^ is the number of ob- 
served jumps, and Cx,...,d is the density of £x,...,d- This result also holds for 
m-dimensional marginal Levy processes with m < d and is stated in Esmaeili 



and Kliippelberg [lOj for two dimensions. 

A straightforward estimation approach would be maximizing the full like- 
lihood function to estimate the dependence structure. This, however, is dis- 
advantageous because of two reasons. The first reason is a numerical one. 
The likelihood function is not easy to evaluate if more than one parameter 
is unknown. The second reason is more conceptual. Since we can only use 
jumps larger than e in all d dimensions, we waste a tremendous part of the 
information about the dependence structure, especially if the dependence 
structure is weak. For weak dependence structures, the probability that two 
jumps are both larger than a threshold (conditioned that at least one jump 
exceeds the threshold) is lower than for strong dependence. 

For both reasons, we estimate the parameters of the bivariate Levy and 
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distributional copulas of the vine structure sequentially. This is also common 
or pair copula constructions of distributional copulas (see, e.g., Hobaek Haff 



12]). We make use of the estimated marginal parameters and start in the first 
tree, using all observations larger than e in the first and second components 
to estimate the parameters of £1^. We continue this procedure for all other 
Levy copulas in the first tree. To estimate the parameter of Ci^, we use 
all observations larger than e in dimensions one, two, and three, as well as 
the previously estimated marginal parameters of the first three dimensions 
and the parameters of (£^2 and £2,3- This means that we proceed tree by 
tree and within the tree, copula by copula or Levy copula by Levy copula, 
respectively. In each step, we make use of the estimated parameters from the 
preceding steps. 

To use the above likelihood for pair Levy copula constructions, we have 
to know how to calculate the density Ci,... )( j of a pair Levy copula. 

Proposition 3. 

Let C\ ... d be a pair Levy copula of the following form 



, Fd\ Z2l ...,z d _ 1 {Ud))dH2,...,d-l{.Z2 } • • • , Zd-l, 

[0,u 2 ]x...x[0,Ud-i] 

and H\...d the corresponding measure and suppose that the density fz...d °f 
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I^2,...,d-i exists. Then the density of exists as well and has the form 



h,...,d{ux, . . . ,u d ) = c(F x 

<9-^l|« 2 ,-,«d-l ^d|U2,...,W d _ 1 

dui dud 

■f2,...,d( U 2i ■ ■ ■ , U d-l)- 

This proposition is proved in the appendix and states that we can iter- 
atively decompose the pair Levy copula into bivariate building blocks and 
therefore evaluate the density function in an efficient manner. 

In contrast to the computation of the density of the pair Levy copula, it 
is not easy to evaluate a higher dimensional pair Levy copula itself. This is 
no real drawback, since the value of £1 d is not needed in most cases. For 
the normalizing constant d of the likelihood, however, £i v .. 5 d has to be 
evaluated. For this step, we apply Monte Carlo methods. The code may be 
obtained from the authors on request, so that for convenience, we omit the 
details here. 

5. Simulation Study 

In order to evaluate the estimators, we conduct a simulation study with a 
5-dimensional PLCC To make the results comparable, all marginal Levy pro- 
cesses are chosen to be stable Levy processes with parameters (a = 0.5, (3 = 
1) and all bivariate Levy copulas in the first tree are Clayton Levy copulas 
(see Equation (j3J)) with parameter 9. The distributional copulas in the higher 
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trees are all Gaussian copulas, i.e., 



where $ p is the distribution function of the bivariate normal distribution 
with correlation parameter p and $ _1 the quantile function of the standard 
normal distribution. 

We analyze three different scenarios of dependence structures: high de- 
pendence (H), medium dependence (M) and low dependence (L). In scenarios 
H and M, we choose a D-vine structure of the PLCC, in scenario L a C-vine 
as it is numerically more appropriate for low dependencies. The D-vine struc- 
ture refers to a structure where all dimensions in the lowest tree form a line 
and are each connected to the nearest neighbors, whereas the dimensions in 
a C-vine structure are connected to only one central dimension (see, e.g., Aas 
et al. [1]). Within a scenario, all Clayton Levy copulas have the same param- 
eter 9 and all Gaussian copulas have the same parameter p. The parameter 
values are summarized in Table [TJ 



Scenario 


Clayton Parameters 6 


Gaussian Parameters p 


High dependence (H) 


5 


0.8 


Medium dependence (M) 


2 


0.3 


Low dependence (L) 


1 


-0.2 



Table 1: Parameters of the PLCC for scenarios H, M and L. 



For each scenario, we simulate a realization of a 5-dimensional Levy pro- 
cess over a time horizon [0,T]. We then estimate the parameters of the 
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process from the simulated data using our estimation approach. We choose 
two different thresholds e = 10~ 4 and e = 10~ 6 for jump sizes we can observe, 
i.e., we neglect jumps smaller than e = 10 -4 or e = 10~ 6 , respectively. Each 
simulation/estimation step is repeated 1000 times. The estimation results 
are reported in Tables [2] and |3j Shown are the true values of the parameters, 
the mean of the estimates of the 1000 repetitions and resulting estimates 
for bias and root mean square error (RMSE). Since the parameters in the 
different trees rely on different numbers of observation (the higher the tree, 
the more dimensions have to exceed the threshold at the same time) we also 
report the mean numbers of available jumps per tree. 

Comparing the two tables, we see that the lower threshold leads to a 
higher number of jumps. We also find that weaker dependence leads to 
less co-jumps available for the estimation of higher trees than a stronger 
dependence. In all cases, the bias is very small. We find, however, that 
the RMSE is affected by the number of jumps available in certain trees as 
it increases with decreasing number of jumps. This effect is illustrated in 
Figure [5] in terms of histograms of the estimates. 

6. Conclusion 

Levy copulas determine the dependence of jumps of Levy processes in a 
multivariate setting with arbitrary numbers of dimensions. In dimensions 
larger than two, however, known parametric Levy copulas are inflexible. In 
this paper, we develop a multidimensional pair construction of Levy copulas 
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Tree # Jumps True Value Mean Bias RMSE 



High Dep. 


1 


870.61 


5 


5.0038 


3.78 


10" 


-3 


2.33 


10" 


-1 




2 


833.51 


0.8 


0.7987 


-1.28 


10" 


-3 


1.33 


10" 






3 


814.39 


0.8 


0.7980 


-1.97 


10" 


-3 


1.34 


10" 


-2 




4 


798.46 


0.8 


0.7890 


-1.10 


10" 


-2 


2.19 


10" 


-2 


Med. Dep. 


1 


707.18 


2 


2.0010 


1.02 


10" 


-3 


9.65 


10" 


-2 




2 


573.56 


0.3 


0.2983 


-1.67 


10" 


-3 


4.58 


10" 


-2 




3 


498.45 


0.3 


0.2983 


-1.72 


io- 


-3 


4.97 


10" 


-2 




4 


451.69 


0.3 


0.3001 


1.31 


10" 


-4 


5.11 


10" 


-2 


Low Dep. 


1 


500.10 


1 


1.0016 


1.63 


10" 


-3 


4.46 


10" 


-2 




2 


267.36 


-0.2 


-0.1987 


1.31 


10" 


-3 


4.98 


10" 


-2 




3 


163.22 


-0.2 


-0.1992 


7.96 


10" 


-4 


7.12 


10' 


-2 




4 


113.91 


-0.2 


-0.2004 


-3.76 


10" 


-4 


9.50 


10" 


-2 



Table 2: Results for a time horizon T=l and a threshold e = 10~ 6 for three scenarios 
from low dependence to high dependence. The columns refer to the number of jumps used 
in the estimation of parameters within a certain tree, the true value of the parameters, 
the mean of the estimated parameters, estimated bias and RMSE from 1000 Monte Carlo 
repetitions. The first three trees contain more than one dependence function and we report 
the mean values of the estimators in these cases. 
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Tree 


# Jumps 


True Value 


Mean 




Bias 


RMSE 


High Dep. 


1 


87.26 


5 


5.0403 


4.03 


io- 2 


7.06 ■ IO" 1 




2 


83.63 


0.8 


0.7933 


-6.76 


10~ 3 


4.61 • 10~ 2 




3 


81.69 


0.8 


0.7810 


-1.90 


io- 2 


5.67- 10~ 2 




4 


80.10 


0.8 


0.7086 


-9.14 


io- 2 


1.46 • HT 1 


Med. Dep. 


1 


70.82 


2 


2.0312 


3.12 


io- 2 


3.19- IO" 1 




2 


57.47 


0.3 


0.2970 


-3.00 


10~ 3 


1.50- IO" 1 




3 


50.00 


0.3 


0.2844 


-1.56 


io- 2 


1.59- IO" 1 




4 


45.37 


0.3 


0.2797 


-2.03 


io- 2 


1.63 • IO" 1 


Low Dep. 


1 


50.21 


1 


1.0246 


2.46 


IO" 2 


1.55 • IO" 1 




2 


26.88 


-0.2 


-0.2019 


-1.87 


io- 3 


1.67- IO" 1 




3 


16.42 


-0.2 


-0.1859 


1.41 


io- 2 


2.57- IO" 1 




4 


11.47 


-0.2 


-0.1378 


6.22 


io- 2 


3.44 ■ IO" 1 



Table 3: Results for a time horizon T=l and a threshold e = 10~ 4 for three scenarios 
from low dependence to high dependence. The columns refer to the number of jumps used 
in the estimation of parameters within a certain tree, the true value of the parameters, 
the mean of the estimated parameters, estimated bias and RMSE from 1000 Monte Carlo 
repetitions. The first three trees contain more than one dependence function and we 
report the mean values of the estimators in these cases. Compared to Tabled the higher 
threshold e results in fewer observed jumps and in higher RMSE of the estimates. 
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High Dependence Setting 



Medium Dependence Setting 



Dependence Setting 




4.E 5 5.5 6 6.5 



1 



0.6 0.7 0.8 0.9 



1 



0.6 0.7 0.8 0.9 



0.6 0.7 0.8 0.9 

e ,5 B 4-° B 




Figure 5: Histograms of the estimation results for a time horizon T=l and a threshold 
e = 10~ 6 . Each column refers to one scenario, the rows refer to the estimated parameters 
in the first to fourth tree. 



from 2-dimensional dependence functions which are either Levy copulas or 
distributional copulas. The resulting parametric Levy copula has the desired 
flexibility, since every regular vine and every bivariate (Levy) copula can be 
used in the PLCC. Applications of the concept can be found in operational 
risk modeling or risk management of insurance companies. In both fields, 
Levy copula models have been proposed but their applicability was limited 
to low-dimensional cases. The pair construction solves these limitations and 
opens the way to high-dimensional applications. In this paper, simulation 
and estimation methods are evaluated in a simulation study. 
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A. Proof of Theorem [2] 

For the proof of Theorem [2J we need a lemma which we state first. 
Lemma 1. 

Let fi be a positive Radon measure on fi : x i— y and f 2 : x i— >■ £| 
fi-measurable measure-valued maps, where and ^ are probability measures 
on R + lozi/i corresponding distribution functions and F^. Let C be a 2- 
dimensional distributional copula and let^ be the probability measure defined 
by the distribution function C(F^,F^) on M^_. Then, the map x >->■ is im- 
measurable. 

Proof: By definition, the maps x h- > ^l(Bi) and x h-> £^(.62) are /x-measurable 
for any Bi,B 2 G B(WL+). This holds in particular for the intervals [0,6] G 
Therefore, the maps x 1— >■ i^(6i) and 2 1— >■ ^(62) are /i-measurable 
for any 6 1; 6 2 G By definition of we have 

tf(S) = C7(F. 1 (6 1 ),F, a (6 a )) 

for any rectangle B G {[0, 61] x [0, 62] | £*i , ^2 G Since C is a copula, 

it is continuous and therefore measurable. We get that x H- £%(B) is a 
composition of /i-measurable functions and therefore /i-measurable for any 
rectangle B G {[0, b]\b G Now that we have shown that x 1— > ^{B) is 

//-measurable for any B G {[0,6]|6 G we use the same argumentation 

as in the proof of Ambrosio, Fusco, and Pallara 2|, Proposition 2.6, to show 
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that x i — y £x{B) i s /i- measurable for any B G B{M? + ). Note that the set 
of intervals B G {[0, b\\b G M^_} is closed under finite intersection, it is a 
generator of the cr-algebra £>(1^_), and there exists a sequence (Bh) of these 
intervals such that M. 2 ^ = UhBh- Denote the family of Borel sets such that 
x i — y is //-measurable by M.. Obviously, D {[0, b]\b G In order 

to use Ambrosio et al. [2], Remark 1.9, we have to show that the following 
conditions hold: 

1. {E h )eM, E h tE^EeM, 

2. E,F, EUFeM^EnFeM, 

3. EeM^ R 2 + \E G M. 

This is already shown in the first part in the proof of Ambrosio, Fusco, and 



Pallara [2], Proposition 2.26. □ 

Now we are able to prove Theorem |2j 

We show that the integral is well-defined in the first step. From Theorem 
[Ufollows that (v,2, . • • , Ud-i) i — y £i|M 2 ,...,u d _i 1S £t2,...,d-i-measurable. By the def- 
inition of measure- valued maps, (m 2 , . . . ,Ud-i) ^i|« 2 ,...,M d _ 1 (-B) is ^2,...,d-v 
measurable for any B G £>(IR + ) and especially for any B G {[0,6]|6 G 
Therefore, 

6|«2,...,Md_i([0) & ]) = - F l|« 2 ,...,n d _ 1 (&) 

is /i2,...,d-i-measurable. With the same arguments, we see immediately that 
-Fd|«2,.--,«d-i(^) i s At2,...,d-i-Hieasurable for any b G Since every copula is 
continuous, we can use the same arguments as in the proof of Lemma [T] to 

34 



show that 



(w 2 , • • • ,tt d _i) !->■ C(Fj 



U2,...,U d . 



-l( u l), F d\ 



U2,...,U d . 



-iW) 



is /i2,...,d-i- m 6asurable and that the integral is well-defined. To show that 
Ci j is indeed a Levy copula, we have to check the properties of Tankov 



22], Definition 3.3. We start by showing that Ci,...^ is d- increasing. In a first 
step, we show this property for any d-box B where all vertices lie in K+. For 
every (uz, . . . , Ud-i) 6 M'jr 2 let £f rf | U2 nrf 1 be the probability measure on M.+ 
defined by the distribution function C(Fi\ U2t ^ jU ._ 1 (wi), Fj,^ ^^(u^)) ■ With 
Lemma[T]we know that («2, . . . , itd-i) H- ^f d | U2 Ud 1 is /i2,... i d-i-measurable. 
By definition of <L\ y .. : d 

£l,...,d(v>l, . . . , Ud) = J C(-^l|2 2 ,. -,zd-i , Fd\z 2 ,...,z d ^ 1 {Ud))dfi2,...,d~l{ z 2, ■ ■ ■ , %d-l) 

[0,M 2 ]X...X[0,M d _l] 



[0,« 2 ]x...x[0,« d _i] \[0,ui]x[0,« d ] / 



d^2,...,d-l{Z2, ■ ■ ■ , z d-l) 



holds, and therefore 



£ w (wi,...,w d ) = ^,...,*-i®^M«a,.-,«d-i([ ' w J x ••• x [ ' n J) 



Mi,...,d([0,wi] x ... x [0,u d ]). 
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Since fi 2 ,...,d-i ® £ 



c 

l,d\u 2 ,...,u d _ 1 



is a positive and well-defined measure, 



V €l _ d (B)= f i 2 _ d . 1 ^^(B)>0. 



In the next step, we denote uj := {u{\i G J} and show that the limit in 
Equation Q) exists for any / C {1, . . . , d} nonempty, I ^ {1, . . . , d}. First, 
suppose that {1, d} C /. Since I ^ {1, . . . , d}, we say w.l.o.g. that {2} ^ I. 
Since €-x,...,d is non-decreasing in every component, it suffices to show that 

lim ...,u d ) 



= lim £i,...,d-i(oo,W2, . ■ ■ ,«d-i) 
u /\{i,d}^°° 

< ^l,...,d-l(00,M2, OO, . . . , OO) = U 2 

to prove that the limes exists. We use the dominated convergence theorem 
(e.g. Ambrosio, Fusco, and Pallara [2|, Theorem 1.21) and the fact that for 
every distributional copula C{u\,u 2 ) < 1 holds. For the inequality, we use 
the fact that Assumption [T] holds for the Levy copula Ci Now, suppose 




Z2,— ,«d_l 



(U d ))dfi 2 ,...,d-l( z 2, ■ ■ -,Zd-l) 




lim 
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that at least one element of {1, d} is not in I. W.l.o.g. {1} ^ I then we have 
lim ...,u d ) 



u i\{d}-* x >J[o )tl2 ]x...x[0,u d - 1 ]J[0,u 1 ] 

= lim / ^.....d-iOza, • • • , 

«l\{d}^oo > /[ 0ini ] x [ 0iU2 ] x ... x [ 0iUd _ 1 ] 

< <£i,...,d-i(«i,oo, ...,oo) = Ui. 

Now that we have shown that the limit exists, it follows immediately that 
£i,.,.,d is (^-increasing on K + . To show that the Levy copula has 
Lebesgue margins, we can again use the same equations as before and re- 
place "<" by "=" since in this case \I\ — d — 1 and therefore we can directly 
use Assumption [TJ □ 

B. Proof of Proposition [2] 

Suppose that i 71 i| U2i ... )?ld _ 1 and ^ are continuously differentiable. 

For any rectangle B — ([0, Ui) x . . . x [0, itj), we get by Theorem [1] 

/ l B {z 1 ,...,z d )dfi K .. td (z 1 ,...,z d )= / Fd\u l ,...,u d ^ 1 (ud)dfii t ... t dr-i(zi, 

JR^_ J[0,u 1 ]x...x[0,u d _ 1 ] 





\z 2 ,...,z d _ 1 {ui), F d \ Z2t ^^ Zd ^ 1 (u d ))dfi 2t ... td -i{z2, • • • , Z d _i 



F 



{ui)dfj, 2 ,..., d ~i(z 2 , . . . , z d -i) 
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By the definition of the pair Levy copula we see that 



/ 1 B (ZI, Zdjdfli^diz!, ...,z d ) 
»/ IR_|_ V J 1R j_ / 

[0,U 2 ]X...X[0,U d _l] 

r 

[0,ita]x...x[0,u d _ 1 ] 

Z2,...,Zd_i 

[0,«i] 



0*11 



z 2 ,...,z d . 



dzi )dfi 2 ,...,d-i( z 2, ■ ■ -,Zd-i) 



[0,u 2 ]x...x[0,ii d _i] 

(*i) 

[0,ui] 



di*i| Z2 ,... A -i(*i) 
5C(Fi| 22i ... i2d _ 1 (21), F dk2i ... i2d _ 1 (u d )) 



[0,ui]x...x[0,u d -i] 



Z 2 ,...,Z d 



22, Zd-iC^l) I d^2,...,d-l( z 2, ■ ■ ■ j 
^A*l,...,(i-l(^l) • • • , ^d-l)i 



and therefore 



5C(-Pl| U 2,...,«d-l( U l)'-^l«2,...,«d-l( M d)) 



9Fi 



u 2 ,...,u d . 



holds /ix j rf_ 1 -almost everywhere. The fact that this result does not only 
hold for fixed values of u d but for all u d G K+ is already shown in the proof 
22] . Lemma 4.2. Since i 71 i| U2 ,...,u d „ 1 , -fd|? i2 ,...,M d _ 1 are continuously 



of Tankov 



differentiable and C is by Assumption [2] also continuously differentiate, we 
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get immediately that Fd\ ul} ... }Ud _ 1 is differentiable and 
dF dl 

ui,...,u d _ 1 { u d) 9 C(Fi\ U2 tUd _ 1 (itl), -fd|u 2 ,...,U d -l M 

is a composition of continuous functions and therefore continuous. Finally, 
all bivariate Levy copulas are by Assumption [2] continuously differentiable 
and therefore, the proposition follows by complete induction. □ 

C. Proof of Proposition [3] 

This statement follows from the definition of the pair Levy copula con- 
struction, since 



(«l) 

[0,« 2 ]x...x[0,it d _i] 



(^l),^,...,^^)) 

[0,m 2 ] x...x [0,«d_i] \ ^ [0,ui] x [0,u d ] 

^-i(^l) ^-x(^) A, f s 

7^ a {Zli z d) J UV>2,...,d-l\Z2i ■ ■ ■ > ^d-lj 

\Z2,-.;Zd-l Z2,--;Zd-l (*0) 



c(-Pl|22,...,* < J_l(2 ; l),-Pd|«J,..., Z d_l(2 ; d)) 

[0,ui]x...x[0,<u d ] 

7^ 77^ J2,...,£i— lV 2 ^? • • • ' z d-l)U{Zl, ■ ■ ■ ,Zd) 



[0,u 2 ] x ...x [0,w d _i] V J [0,«i] x [0,u d ] 
/2,...,d-l(^2, • • • , ^d-l)d(^2, • • • , Zd-l) 
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as stated. □ 



40 



References 

[1] Kjersti Aas, Claudia Czado, Arnoldo Frigessi, and Henrik Bakken. Pair- 
copula constructions of multiple dependence. Insurance: Mathematics 
and Economics, 44(2):182-198, 2009. 

[2] Luigi Ambrosio, Nicola Fusco, and Diego Pallara. Functions of Bounded 
Variation and Free Discontinuity Problems. Oxford University Press, 
2000. 

[3] I.V. Basawa and P.J. Brockwell. Inference for gamma and stable pro- 
cesses. Biometnka, 65(1):129-133, 1978. 

[4] I.V. Basawa and P.J. Brockwell. A note on estimation for gamma and 
stable processes. Biometrika, 67(l):234-236, 1980. 

[5] Tim Bedford and Roger M. Cooke. Probability density decomposition 
for conditionally dependent random variables modeled by vines. Annals 
of Mathematics and Artificial Intelligence, 32(l-4):245-268, 2001. 

[6] Tim Bedford and Roger M. Cooke. Vines: A new graphical model for 
dependent random variables. The Annals of Statistics, 30(4):1031-1068, 
2002. 

[7] Daniel Berg and Kjersti Aas. Models for construction of multivariate 
dependence: A comparison study. The European Journal of Finance, 
15(7):639-659, 2009. 



41 



[8] R. Cont and P. Tankov. Financial modelling with jump processes. 2004. 

[9] Habib Esmaeili and Claudia Kliippelberg. Parameter estimation of a 
bivariate compound poisson process. Insurance: Mathematics and Eco- 
nomics, 47(2):224-233, 2010. 

[10] Habib Esmaeili and Claudia Kliippelberg. Parametric estimation of 
a bivariate stable Levy process. Journal of Multivariate Analysis, 
102(5):918-930, 2011. 

[11] Habib Esmaeili and Claudia Kliippelberg. Two-step estimation of a 
multivariate Levy process, preprint, 2011. 

[12] Ingrid Hobaek Haff. Comparison of estimators for pair-copula construc- 
tions. Journal of Multivariate Analysis, forthcoming, 2012. 

[13] Harry Joe. Families of m-variate distributions with given margins and 
m(m — l)/2 bivariate dependence parameters. Distributions with Fixed 
Marginals and Related Topics, pages 120-141, 1996. 

[14] M. Johannes. The economic and statistical role of jumps in continuous- 
time interest rate models. The Journal of Finance, 59:227-260, 2004. 

[15] Olav Kallenberg. Foundation of Modern Probability. Springer, 2002. 

[16] Jan Kallsen and Peter Tankov. Characterization of dependence of multi- 
dimensional Levy processes using Levy copulas. Journal of Multivariate 
Analysis, 97:1551-1572, 2006. 



[17] John Kingman and Samuel J. Taylor. Introduction to measure and prob- 
ability. Cambridge University Press, Cambridge, 1966. 

[18] Suzanne S. Lee and Jan Hannig. Detecting jumps from Levy jump dif- 
fusion processes. Journal of Financial Economics, 96(2):271-290, 2010. 

[19] R.B. Nelsen. An Introduction to Copulas. Springer, New York, 2nd 
edition, 2006. 

[20] J. Rosinski. Levy Processes - Theory and Applications. In: Series Rep- 
resentation of Levy Processes from the Perspective of Point Processes. 
Birkhauser, Boston, 2001. 

[21] Ken-Iti Sato. Levy Processes and Infinitely Divisible Distributions. Cam- 
bridge University Press, Cambridge, 1999. 

[22] Peter Tankov. Simulation and option pricing in Levy copula models. 
Unpublished manuscript, 2005. 



43 



